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ABSTRACT 

We analyse the BATSE 3B catalog using the pair-matching statistic. This 
statistic counts only the burst pairs which may have originated from the same 
source, so it is less likely to yield false detections of "repeating bursts" than the 
nearest neighbor and correlation function statistics. Even in the ideal case when 
repeating is the only possible source of burst correlations, the pair matching 
statistic is more sensitive to repeating bursts than these other statistics particu- 
larly for models which predict faint or multiple burst repetitions. We find that 
the BATSE 3B data set contains no excess of matched burst pairs over the expec- 
tation from a sample with random positions. We also apply the pair-matching 
statistic to the bursts that previously appeared in the BATSE IB catalog which 
now have improved positions and position errors in the BATSE 3B data set. 
Previously, these bursts had exhibited some peculiar position correlations that 
were interpreted by some as evidence for burst repetition, but we find that these 
correlations have disappeared with the improved BATSE 3B positions. 



Submitted to The Astrophysical Journal Letters 



1. Introduction 



Whether classical gamma ray bm"Sts might recur from the same somxes on 
the sky is an important issue because it has direct bearing on the nature of the 
burst sources: Do the gamma ray bursts originate from from inside the Galaxy 
or from cosmological distances? In cosmological models, the bursts appear to be 
distributed uniformly on the sky and the observer appears to be at the center 
of a distribution which has a dearth of faint (and presumably distant) bursts 
due to cosmological redshift elfects. These predictions are consistent with the 
burst distribution observed by BATSE (Meegan, et ai, 1992; Fishman, et ai, 
1993, Meegan, et ai, 1995a). On the other hand, the discovery of repeating 
bursts would be an evidence against cosmological models because they generally 
demand an energy output so large that the source must be consumed during the 
burst. Clearly, such a source can not repeat. 

The implications of the publicly released BATSE IB data set regarding burst 
repetition have been controversial. Quashnock and Lamb (1993) presented a 
"nearest neighbor" analysis of this data and claimed statistical evidence for 
repeating bursters. Other authors (Narayan and Piran, 1993; Nowak 1994; 
Strohmayer, Fenimore, and Miralles, 1994; Meegan, et al., 1995b; and Efron 
and Petrosian 1995) have also analyzed the BATSE data primarily with nearest 
neighbor statistic or the angular correlation function. The correlation function 
analyses have found an excess of pairs at separations close to 0° as well as a sim- 
ilar excess at 180°. The excess at 180° is difficult to explain as a physical efi^ect, 
so perhaps the correlations at both 0° and 180° are due to a systematic error in 
the burst locations. One peculiarity noticed by many of these authors was that 
the "signal" in the correlation function occurred on a scale of ^ 4° while the 
estimated position errors of the bursts responsible for the correlations were ^ 8°. 
Quashnock and Lamb suggested that this was the signal of a burst population 
with a few repeating bursts that were seen several times. It was also noted that 
many of the bursts responsible for this "signal" were relatively faint. 

One problem that affects both the nearest neighbor statistic and the two-point 
angular correlation function is that these statistics really measure correlations 
between burst positions rather than the property that we are really after: whether 
each pair of bursts could have come from the same source. It is true that repeating 
burst sources will cause burst positions to be correlated, but there are other 
mechanisms that might also cause the bursts to be correlated, such as systematic 
position measurement errors. We need a statistic that will not only measure burst 
correlations but will determine whether each pair of bursts could have come from 
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the same source. If there is a significant excess of 'matched' burst pairs, then 
is the separation distribution consistent with the expected distribution from a 
population of repeating sources? 

In order to address this problem, we have devised the pair matching statistic^ 
(Bennett & Rhie, 1996: Rhie & Bennett, 1995). For each pair of bursts, we de- 
fine a match probability that measures the likelihood that this pair has originated 
from the same location on the sky. The match probabilities are then combined 
into our pair matching statistic which measures only those burst pairs which are 
consistent with being produced by the same source. Even in the idealized case in 
which all burst correlations are due to burst repetition, the pair matching statis- 
tic is more sensitive to burst repetition than the nearest neighbor and correlation 
function statistics particularly for faint or multiple repeating bursts. The sensi- 
tivity of the pair matching statistic to the fraction of repeating pairs improves 
as 1/N where TV is total number of bursts. (The noise scales as the square root 
of the total number of burst pairs). In contrast, the sensitivity of the nearest 
neighbor statistic improves much more slowly than this for large N. 

Another important distinction is due to the fact that the pair-matching statis- 
tic makes use of the match probability for each pair of bursts. This allows us 
to compare the shape of the observed pair match probability distribution to the 
distribution predicted by repeater models. Thus, if the pair matching statistic 
for the BATSE data should reject the "null" hypothesis of no repeating bursters, 
then we can also apply a "goodness of fit" test to see if the observed signal is 
actually consistent with burst repetition. In fact, when we applied (Bennett & 
Rhie, 1996) the pair-matching statistics to the BATSE IB and 2B data sets, we 
found that the observed small angle correlation was due to bursts separated by an 
angle smaller than expected based on the BATSE position error estimates. When 
all bursts pairs consistent with originating from the same source were considered, 
no significant signal of repeating bursts remained. 

The main point of this paper is to apply our pair-matching statistic to the 
bursts in the BATSE 3B catalog. In section 2, we define the pair matching 
statistic and summarize our previous results. In section 3, we analyze the BATSE 
3B with the pair matching statistic, and we also consider the BATSE IB' data 
set which consists of the BATSE IB bursts with improved BATSE 3B positions. 



t A similar statistic has been developed by Petrosian & Efron (1995) to study repeating on 
short timescales. 
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2. The Pair Matching Statistic 



The pah matching statistic is discussed at length and compared to other 
statistics in our previous paper (Bennett & Rhie, 1996), so we will just summarize 
it here. If we let 6ij be the angular distance between bursts i and j, then the joint 
probability density of the pair is e~^ij^'^"ij where afj = crfg^g^^ + crfsys + 
^jstat + ^jsys quadrature sum of the statistical and systematic position 

errors for bursts i and j. The match probability for bursts i and j is defined to 
be 

Pij = 4" / e"^'/^''^-'- ed0 = e-^?'/^^- . (2.1) 

which is just the formal probability that if bursts i and j originated from the 
same position on the sky, they would have a separation at least as large as that 
observed. Thus, for burst pairs that come from the same source, we should expect 
Pij to be uniformly distributed between and 1. 

One can express the expected distribution of pij values as a distribution 
function: g{p). For matched pairs, g{p) is just a constant, so for iVmatched matched 
pairs, we have 

^matched (P) = -^matched , (2-2) 

while for a set of -A^random randomly distributed burst pairs we have 

„ erf - ^random (cr^^y•, 

^random (P) = 2^ 3 ^ 2 ' ^^'"^^ 

i,j<i 

where we have used the small angle approximation to derive eq. (2.3). Note 
that (/random (p) does not actually diverge for small p because there is a minimum 
allowed value of p which corresponds to = tt. In order to test for repeating 
bursts, we will need to be able to detect a signal of the form, (2.2), from the 
background noise of the form (2.3) with iVj-andom ^ -^matched- Using eq. (2.1), 
we can define the pair match statistic: 

Qa{p) = Y.P'ij^^Pi3-p) ^ (2-4) 
pairs 

where G is the familiar step function. The value of Qaip) for an uncorrelated 
data set depends on the number of data points and the distribution of error bars. 
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so Monte Carlo simulations of random data sets with the observed distribution 
of error bars are required to assess the statistical significance of these statistics. 
For a gamma ray burst population that does repeat, then the expected signal is 
given by 

Qaip) - (Qa(p))KANDOM = ^^ff (l "P'^^^) " E 4 ' (2-5) 

matched 
pairs 

where (<3q(p)) random mean value of the Monte Carlo simulated BATSE 

catalogs with random positions, and A^matched is the number of burst pairs origi- 
nating from the same physical burst location on the sky. For a = 0, (1 — p°')/a 
is replaced by — Inp in the last term of eq. (2.5). In our previous paper, we 
compared the expected signal-to-noise ratio for all values of p and a = or 1. In 
all cases, we found that Qi{0) was the most sensitive statistic, so this the statistic 
that we use determine if a repeater signal can be seen in the data. 

We have also compared the sensitivity of the pair matching statistic to that of 
the two-point correlation function and the nearest neighbor statistics and found 
that the pair matching statistic is generally more sensitive to repeating bursters. 
In contrast, the nearest neighbor statistic is particularly insensitive to repeating 
in large data samples: In a sample of 1122 bursts with 240 burst pairs from 
the same location, the nearest neighbor statistic will only see a signal which is 
significant at the 60-90% confidence level while the pair matching statistic sees 
a signal which is significant at the 3-cr level (assuming repeater burst position 
errors of 8°). 

We should point out that although the value of p in eqns. (2.1) and (2.4) is 
a calculated assuming Gaussian errors, the pair matching statistic does not rely 
upon the assumption that the errors are Gaussian distributed. To the extent 
that the position errors do not follow a Gaussian distribution, we can consider 
Pij to be a parameter which is correlated with rather than equal to the true 
match probability. Since we are comparing against the "null hypothesis" of no 
repeating burst sources, our Monte Carlo simulated BATSE catalogs use only 
random positions. Thus, the comparison to the random catalogs remains valid 
for any distribution of measured position errors. 
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3. Pair Matching Analysis of the BATSE 3B Dataset 



Before we present the pair matching results for the BATSE 3B data set, let 
us discuss the subsets of the data that we will apply the statistic to. Previous 
authors (Quashnock and Lamb, 1993; Narayan and Piran, 1993) have considered 
sub-samples of the catalog defined by the requirement that crgtat be less than 
some value (9° for QL and 4° for NP). For comparison with their results and 
with our BATSE 2B results, we also include these sub-samples in the present 
analysis. Out of 1122 bursts in the BATSE 3B catalog, 956 and 628 bursts pass 
the (Jstat < 9° and the crgtat < 4° cuts respectively. The BATSE IB' data set is 
the subset of the BATSE 3B data set which corresponds to the same time interval 
as the BATSE IB data set. Of the 263 bursts in this data set, 229 bursts pass 
the cTgtat < 9° cut and 147 bursts pass the dstat < 4° cut. In the original analysis 
of the 260 BATSE IB bursts, 202 and 133 bursts passed these two cuts. These 
differences are due to burst reclassification and improved positions. 

With data sets containing a few hundred to a thousand bursts with position 
errors ranging from a few to a few tens of degrees, accidental matches between 
bursts from different sources will be common, and they will be a source of noise 
for the pair matching statistic. The number of accidental matches will depend 
somewhat on the average exposure of BATSE to different areas of sky. For exam- 
ple, if BATSE observed only the Northern hemisphere, the number of accidental 
matches would be increased because the observed bursts would be confined to a 
smaller region of sky. We correct for BATSE's non-uniform exposure by using 
the the exposure table provided with the BATSE data. However, these exposure 
tables do not completely describe the exposure distribution for the faintest bursts 
that BATSE can detect. For this reason, we have defined on exposure corrected 
sample which contains only those bursts for which the sky exposure can be accu- 
rately determined from the the BATSE exposure table. The details of this event 
selection are given in Bennett and Rhie (1996). There are 540 BATSE 3B bursts 
which pass this cut. 151 of these are in the BATSE IB' data set compared to 
150 in the BATSE IB data set. 

Fig. 1 shows the variation of the Qi{p) statistic with p for data samples 
from the BATSE 3B catalog. We have chosen to plot the differential version of 
this statistic {i.e.AQi(p) = Qi{p) — Qi{p + Ap)) because this makes it easier 
to see the p dependence of Qi{p). Also plotted are the expected signal from a 
sample which has 4 x 10~^ of all burst pairs originating from the same source 
and the Monte Carlo average of data sets with random positions. Clearly, only 
the exposure corrected sample might have a signal as large as ~ 4 x 10~^ level. 
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and as shown in table 1, number of matched pairs for this sample is only about 
1 a above the average of the Monte Carlo simulations. 

The precise limits on the number of burst pairs that might originate from 
the same source are shown in Table 1. For the full sample the 2 o upper limit on 
the number of bursts that might have originated from the same sources is 120 or 
about 1.9 X 10~^ of all the burst pairs. It is somewhat more conventional to quote 
the number of bursts that are due to repeating sources, but this is an ambiguous 
number for the pair matching statistic (as well as for the correlation function) 
because it is based on burst pairs rather than single bursts. For example, if each 
burst is seen to repeat only once, a fraction of 1.9 x 10"'^ repeat pairs corresponds 
to 240 single bursts or 21 % of the 1122 total that are from repeating sources. 
However, the 120 pairs could also be due to 1 source that generates 16 bursts. 
In this case, the fraction of all bursts due to repeaters would be 1.4%. 

Since our main question is whether there is any statistically significant sig- 
nal of burst recurrence in the same position, the fraction of simulated BATSE 
catalogs with random burst positions that have Qi(p) > the observed value is of 
particular interest. These values are plotted in Fig. 2 for both the BATSE 3B 
data set and its IB' subset for both "matches" and "antipode matches'.' These 
plots indicate that there is no significant signal for any of our data subsets 

for any value of p except for perhaps the BATSE 3B antipodal pairs for p > 0.7 
where the random probability of exceeding the observed value is ^ 0.01. 

Given the large number of data subsets and the range of p values we've examined, 
this can hardly be considered significant, however. 

It is of particular interest to note that the BATSE IB' data set does not show 
any hint of a matched pair or antipode pair signal. This differs from the results 
for the original release of the BATSE IB data set. In the original IB data set, the 
"match" statistics for the cTgtat < 9° sample (and to a lesser extent in some of the 
other samples), there appeared to be some sort of a signal in the measured Qi(p) 
values for p ^ 0.6. The random probability of the observed values was a few 
xlO~^, but for smaller p values where the signal to noise should have been larger 
the random probability of the observed value rose to reach about 2% at p = 0. 
Thus, there appeared to be a significant excess of burst pairs separated by angles 
smaller than expected for repeating bursts and a much more marginal signal 
that might have been consistent with a small population of repeating bursts. 
Undoubtedly, burst pair excess at p ^ 0.6 is the same signal that had been seen 
by other authors using the nearest neighbor and correlation function statistics. 
Thus, the refined burst position analysis used for the BATSE 3B data release has 
clearly removed the position correlations that some had interpreted as a signal of 
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repeating bursts. The peculiar correlation at 180 degrees has also disappeared. 
This suggests that both the small angle and the antipodal clustering seen in the 
original BATSE IB data set was due to systematic effects associated with the 
method of determining burst positions. 

In summary, we have used our pair matching statistic to search for evidence 
of burst repetition in the BATSE 3B data set, and we find no evidence for burst 
repetition at a level of sensitivity that is superior to the sensitivity that can 
be attained with the nearest neighbor statistic or the correlation function. We 
have also analyzed the BATSE IB' subset of the BATSE 3B data set and have 
found that the position correlations that existed in the IB data set have now 
disappeared with the improved BATSE 3B position measurements. This suggests 
that the correlations in the BATSE IB data might have been caused by systematic 
errors in the original position determinations that have now been corrected. 
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Sample: 


All 


""Stat < 9° 


(Tstat < 4° 


Exposure corrected 


-^match 3B 


-43 ± 82 


-85 ± 47 


-17± 20 


32 ±29 


-^antipode 3B 


25 ±82 


59 ±47 


19 ±20 


-18 ±29 


match significance 3B 


68.55% 


96.74% 


81.50% 


13.74% 


antipode significance 3B 


36.40% 


10.51% 


16.35% 


72.35% 


-^match IB' 


2± 19 


5± 12 


2±5 


4±7 


-^antipode IB' 


-8± 19 


10± 12 


4±5 


9±7 


match significance IB' 


45.15% 


33.61% 


30.52% 


71.45% 


antipode significance IB' 


65.57% 


19.17% 


20.83% 


9.51% 


-^match IB 


13± 25 


27± 14 


1 ±7 


3± 12 


-^antipode IB 


54± 25 


24± 14 


17 ±7 


9± 12 


match significance IB 


29.3 % 


2.65% 


43.36% 


36.96% 


antipode significance IB 


2.42% 


4.54% 


1.32% 


21.92% 



Table 1. The number of matches and anti-matches is given for the BATSE 3B, 
IB', and IB data sets as measured with the Qi(0) statistic with 1-cr error bars 
determined from the RMS deviation from the mean of 10,000 simulated BATSE 
catalogs. The significance refers to the fraction of the 10,000 simulated catalogs 
with larger values of (5i(0) than the real data. The IB' data set is the subset of 
the 3B catalog that corresponds to the IB catalog. 
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FIGURE CAPTIONS 



1. AQi{p) = Qi{p) - Qi{p + Ap) for the BATSE 3B catalog is plotted as a 
function of p for all bursts, (a), bursts with o"stat < 9°, (b), bursts with 
Cstat < 4°, (c), and bursts which pass the "exposure correction" cut, (d). 
Ap = 0.1 is the histogram bin size. Qi{p) is obtained from AQi{p') by 
summing over all p' > p. 

2. The fraction of 10,000 simulated BATSE datasets which have Qi{p) values 
which exceed Qi{p) for the real data is plotted as a function of p for (a) 
the BATSE IB' subset of the BATSE 3B catalog, and (b) the BATSE 3B 
catalog. 
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